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The fate of R-parity is one of the central issues in the minimal supersymmetric standard model 
(MSSM). Gauged B — L symmetry provides a natural framework for addressing this question. 
Recently, it was pointed out that the minimal such theory does not need any additional Higgs if 
the B — L breaking is achieved through the VEVs of right-handed sneutrinos, which ties the new 
physics scale to the scale of the MSSM. We show here that this immediately leads to an important 
prediction of two light sterile neutrinos, which can play a significant role in the BBN and neutrino 
oscillations. We also discuss some new relevant phenomenology for the LHC, in the context of the 
minimal supersymmetric left-right symmetric theory which provides a natural setting for the gauged 
B ~ L symmetry. 



I. INTRODUCTION 

Over the decades, low energy supersymmetry became 
the central extension of the standard model (SM). Its 
attributes are well known: stability of the gauge hierar- 
chy, gauge coupling unification [l| and the radiative Higgs 
mechanism However, supersymmetry allows for new 
interactions which break the B and L accidental sym- 
metries of the SM. The stability of the proton and the 
associated phenomena become a mystery. Due to non- 
renormalization theorems, one can simply make the new 
interactions (or a subset of them) vanish. Only if all of 
them end up vanishing, one gets the well-known R-parity 
symmetry, under which ordinary particles are even and 
super-partners are odd. R-parity is equivalent to matter 
parity 

A/p = (-1)3(5-^) . (1) 

A priori, there is no way of knowing whether or not R- 
parity is a good symmetry. Gauging B — L symmetry 
sheds light [3| on the fate of this otherwise ad hoc sym- 
metry. This can be achieved at the expense of introduc- 
ing right-handed neutrinos, which is the simplest way 
of generating neutrinos masses. The right-handed neu- 
trinos are particularly natural in the context of left-right 
(LR) symmetric theories Q , which provide a possible un- 
derstanding for parity violation in nature. Furthermore, 
they contain automatically gauged B — L symmetry. In 
short, these theories offer a unified picture of three fun- 
damental issues: the origin of parity violation, neutrino 
mass and the fate of R-parity. 

Recently new progress was made in the construction 
of these theories, by achieving the breakdown oi B — L 
and/or LR symmetries through the RH sneutrino VEVs 
only 0, @ , which simultaneously breaks R-parity. These 
are arguably minimal such models and furthermore, they 
tie automatically the gauge symmetry breaking scale to 
the scale of supersymmetry breaking soft terms. This 
provides a great boost towards the detectability of this 
class of theories at the large hadron collider (LHC). 



We show here, that these theories have an additional 
important feature: the necessary existence of two light 
(sub-eV) sterile neutrinos which may play a crucial role 
in the anomalies of neutrino oscillations. Furthermore, 
due to the new gauge interactions, they are equally abun- 
dant as the left-handed ones at the time of big bang nu- 
cleosynthesis (BBN). This is the central finding of this 
Letter. The latest analysis of BBN seems to favor such 
degrees of freedom Q. We also discuss some new asso- 
ciated phenomenology for the LHC in the context of the 
LR symmetric version of the theory. 

Before turning to the description of the model, let us 
clarify the essence of our result. First of all, what role 
does supersymmetry play? After all, gauging B — L is 
interesting in its own right and it is worth discussing 
what happens without supersymmetry. Obviously, the 
neutrino spectrum depends crucially on what Higgs one 
chooses in order to break B — L. In the seesaw picture [1], 
one opts for the B — L = 2 SM singlet field. One may 
instead prefer the B — L = 1 Higgs field, in which case 
at the renormalizable level, the RH neutrinos do not get 
Majorana masses and one is left with 3 Dirac neutrinos. 
This was discussed in Ref. Q. The modern point of view, 
though, is not to ignore higher dimensional operators, as 
illustrated by the famous d = 5 Weinberg operator in the 
SM [l^ . The analogous operator in this case implies the 
RH neutrino masses on the order of Mgj^/M, where M 
is the scale of new physics above Mbl- Both of these 
scales are arbitrary and no prediction can be made. One 
could also choose a large B — L charge (> 3) for the SM 
singlet Higgs, in which case, RH neutrinos would remain 
massless and again one would end up with three Dirac 
neutrinos. In short, the resulting physics depends on 
the arbitrary model building and/or assumptions about 
higher dimensional operators, with no information about 
mass scales. 

On the other hand, as we argued, the supersymmetric 
B — L theory in its minimal form relates the scale oi B~L 
breaking to the scale of low energy supersymmetry, and 
leads automatically to two extra sterile neutrinos. In 



what follows, we show now how this takes place. 



II. THE SUPERSYMMETRIC [/(I)' MODEL 

The minimal model [Hj of gauged B — L symmetry is 
based on the gauge group ^ Gb-l = SU{3)c x SU{2)l x 
U{1)y X U{1)b-l- The matter assignment is given by 



Q 



L 



(2) 



where beyond the usual SM particles, one needs (anti) 
RH neutrinos, one per generation, in order to cancel the 
B — L anomaly. The most general gauge-invariant super- 
potential becomes then (we ignore the quark sector being 
irrelevant for what follows) 



(3) 



where Hu,Hd are the MSSM Higgs doublets. In this 
minimal model, one does not introduce any new Higgs 
multiplets so that B — L symmetry must be broken by 
the non- vanishing VEV of the RH sneutrino field {N'^). 
This requires a tachyonic soft mass for the N'^ field. 

The most appealing feature of this approach is that the 
scale oi B — L breaking is tied to the scale of soft super- 
symmetric terms, and thus if the MSSM is to be observed 
at LHC, so would be the new gauge interactions. An- 
other important feature is the prediction of light sterile 
neutrinos, as we show now. 



A. Neutrino Masses at Tree Level 

The essential point is one RH neutrino becomes mas- 
sive at the first stage of symmetry breaking. One is al- 
ways free to rotate the RH sneutrino VEVs into the di- 
rection (Ni). In turn one gets 



M 



N?,Zbl 



Mz' 

MZ' TOl/2 



A/^c = A/«c = 0, (4) 



where Mz' = 9bl{Ni) is the mass of the Z' boson, Zbl 
is the gaugino of U{1)b-l with generic Majorana soft 
mass mi/ 2 and gBL is the B — L gauge coupling. For 
simplicity, from now onwards, we will assume roughly 
one scale of soft supersymmetry breaking terms, which is 
denoted as Msusy, thus mi/2 — Msusy ■ 

After the second stage of symmetry breaking where 
Hu,Hd gets their VEVs, from Eq. (jS]), one gets Dirac 
neutrino mass terms. One linear combination of LH neu- 
trinos, called vi hereafter, get a Majorana mass m^-^ ~ 



{Yj-^)li{Hu) /M]sr^ through the seesaw mechanism, when 
Ni gets integrated out ^. To implement TeV seesaw, we 
need {Yj^)ii < 10~^. In principle, the left-handed sneu- 
trino ly can also get a VEV, which has to be very small. 
The point is that it mixes LH neutrinos with neutralinos. 
A different linear combination of LH neutrinos, called 1/2 
hereafter, gets a Majorana mass to^^ oc {v)"^ / 1711/2, when 
neutralinos are integrated out. Neutrino oscillations to- 
gether with beta decay require the mass eigenstates to 
lie below eV. For mi/2 ^ 100 GeV - ITeV, from < 
eV, one gets (v) < 300 KeV- IMeV. 

In general, the five light neutrino states 



1^2,1^3, Ni,NI) 



will mix and produce the follow- 



ing mass matrix at renormalizable level 
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where the state 1^3 is orthogonal to i>i , 1^2 states and ~ 
Y^ (Hu) are the Dirac masses in this basis. At tree-level, 
there are five Majorana light states. Their masses are not 
predicted and must be inferred from experiment. There 
is the cosmological bound from the WMAP on the sum of 
all light neutrino masses, again on the order of eV. This 
implies the corresponding {Yj^)ij < 10^^^ for i = 1,2,3 
and j = 2, 3. 

The extra two light states may be welcome [ll| . in 
view of the oscillation anomalies in LSND, MiniBooNE 
and MINOS (l^ - [l^ . Furthermore, the case is being 
made precisely for two more light fermionic states, in 
order to have successful BBN Q- These states by defi- 
nition, should be roughly equally abundant as the usual 
three active neutrino in the SVI. This is guaranteed due 
to the B — L gauge interactions, which keep the light 
states N2,N^ in thermal equilibrium until the temper- 
ature falls down to about lOMeV. In principle Yukawa 
coupling could also do the job |15| . which would require 
a quantitative analysis, not needed here. 



B. Higher Dimensional Operators 

How robust are these predictions? What happens if 
one introduces higher dimensional operators? The gauge 
singlet combinations which are building blocks for such 
operators are 



LHuN', HuHd, LH^N', N'^N' 



(6) 



^ Strictly speaking, one can start with two U (1) symmetries, which 
are arbitrary linear combinations of B — L and hypercharge. The 
difference lies in the couplings of Z' boson to the SM fermions. 



^ Notice there is another contribution due to the LH neutrino- 
Higgsino-neutral gaugino mixing, which gives Majorana mass to 
the same linear combination v\. 
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The first two appear in the supcrpotcntial at tree level 
and thus offer nothing new. In this sense they can be clas- 
sified as trivial operators, since they are basically rcnor- 
malizing the already present couplings and masses. The 
typical operators that we have in mind can be represented 
by the following F-type term 



Of = [LHuN'^) 



HuHd 



(7) 



This would induce a Dirac mass term of the order of 
mo — M^/AP. Clearly, cither the scale of new physics 
is large enough, M > 3 x 10^ GeV or it requires a small 
Wilson coefficient. 

Similarly, the third term in Eq. (6) also gives a Dirac 
mass term from Kahler potential, mo = {Fh^)"^ /AI = 
li{Hu)/M ~ MwMblIMv[\ which now requires A/ > 
10^'* GcV. This is in complete analogy with the situation 
in the SM, where d = 5 Weinberg operator requires a 
large scale or a small coupling. 

Notice this is nothing new. In the SM, the tree level 
Yukawa couplings are similarly affected by {W H / A'P)" 
Wilson-like terms. Small Yukawa couplings as above im- 
ply large M or small coefficients. This is the cornerstone 
of the SM flavor structure. For this reason we do not 
discuss such operators anymore. 

More interesting is the fourth type gauge invariant op- 
erator in Eq (6), which affects the tree- level Kahler po- 
tential. Together with HuHd, it produces a general class 
of operators relevant for RH neutrino Majorana masses ^ 



O 



(p.<?) 

D 



M2 



HuHd 



'hIhI^ 



(8) 



where a,b ~ 2,3. The case p = q = leads to the sterile 
neutrino mass matrix (corresponding to the lower-right 
2x2 submatrix of Eq. (5)) 



m^l(0,0) 



M2 



Y, 



D 



M2 



< lO^^eV, 
(9) 

for M > lOMz'- The reason this comes out negligible is 
the double suppression of small Yj^ and small (57) . 

Ironically, the leading contribution comes from the 
next order operator with p = 0, q ^ 1, which now leads 
to 

„t(0, 1) ^ '"■'".^■:>«' . < MOV , 



A/4 



Af4 



(10) 

for M > lOMz'- It is much larger than m^(0, 0) since 

If these sterile neutrinos N2.3 are to have appreciable 
masses, they must decay. The lighter of the two (called 



^ In principle, the operators contain additional factors 
{HlHu)"^{HlHa)'^, but they clearly lead to same physics. 



Nl hereafter) decays through the dipole operator, which 
is generated at the loop level. 



(11) 



while the heavier one (called Nh) can also decay via the 
higher dimensional operator from the Kahler 



L^LNjNH 



AP 



(12) 



It is easy to verify that both lifetimes turn to be longer 
than 10^2 sec, for m^v^ <MeV, or cosmologically sta- 
ble for mj^^ < 100 KeV. Since A^^, Nh are roughly as 
abundant as the photons, unless being very light, they 
would quickly dominate the universe when the photon 
temperature drops below their mass. This would run 
into conflict with observations on BBN, CMB and large- 
scale structure. Recent studies of the constraints from 
structure formation find to^ ^ 1 eV [l^ , and in turn 
A/ > 3 X WMw 

What about the loops? They arc actually very tiny if 
nonzero, for they must go through the Dirac couplings 
Y^ , which is the only source of the breaking of flavor 
symmetry in the neutrino sector. 

In summary, should the future confirm the need for 
sterile neutrino masses in the eV or sub-eV region, the 
scale Af of new physics beyond B — L would have to be 
fairly low, on the order of 10'"' GcV. This is good news 
since any oasis in the desert is welcome. 



III. LEFT-RIGHT SYMMETRIC VERSION 

The minimal supcrsymmctric LR model (which we call 
the FS model Q here and below) is based on the gauge 
group Glr = Si7(3)e x SU{2)l x SU{2)r x U{1)b-l- Be- 
sides the Q and L doublets, one has the singlet fermions 
of Eq. (O as the SU{2)r doublets 



0^ = 



(13) 



Now Hu and Hd form a bidoublet <i> under SU[2)l x 
SU{2)ii and as before no new Higgs field is introduced. 
Again R-parity is a gauge symmetry and is spontaneously 
broken by the VEV of RH sneutrino. This requires the 
asymmetric soft terms for left and right sleptons, as as- 
sumed in 0. Strictly speaking, this is not a left-right 
symmetric theory although the breaking is soft. In prin- 
ciple, this soft breaking can emerge from the D-term VEV 
of a parity odd singlet (partially) responsible for super- 
symmetry breaking in the hidden sector. 

The characteristic of the model is the prediction for 
the ratio of heavy neutral and charge gauge bosons 



A//z' 



Af, 



Wa ^JigR/ghY - tan2 Ow 



> 1 



(14) 
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where 0w is the Weinberg angle. Besides the usual upper 
limit on gpi from the pcrturbativity requirement, it is also 
easy to see there is a lower limit, Qr > gh tan^^/- When 
9l ~ gB., one gets Mz'/Mw^ = cos 6w / y/cos 29w ~ 1-2. 
It should be contrasted with the popular version of the 
theory based on the rcnormalizable seesaw mechanism. 
In the latter case, one utilizes the triplet Higgs to break 
the SU{2)[i symmetry and give mass to all TV^, and one 
gets Mz'/Mwn w 1.7 for ql « qr. 

A. Neutrino Masses 



FIG. 1: The Cheops diagram: potentially large radiative cor- 
rection to the sterile neutrino masses, where a,b — 2, 3. 



What about the neutrino mass spectrum in this the- 
ory? In the same manner as in the U{lY model, one ends 
up with only one heavy RH neutrino at tree level. The 
point is that we only added RH charged current interac- 
tions without modifying the neutral fermion mass matrix 
structure. Regarding the higher dimensional operators, 
the conclusion goes through as before, which again im- 
plies two light sterile neutrinos. The only difference lies 
in the loop corrections to the sterile neutrino masses. 
There is a new contribution shown in Fig. 1. 

The point is that in the LR model, there are bidoublets, 
and so the i/r coupling might enter the above loops. The 
most general superpotential takes the form 

W = YeL<PiL^ + y,^L$2^" + MuTr[$,$,] , (15) 

where — {Hu,Hd)i, i,j = 1,2, with the most general 
VEVs 

In this case, the Dirac neutrino mass is still proportional 
to y,^, with (2/^)n ~ 10-6, < 10-", = 2,3), 

similar to the B-L model. The main contribution to 
charged lepton mass is from the first Yukawa Yi = 
diag(ye, JZ/j, J/t) and the second term is only a small cor- 
rection. A crucial point to note is that Hui has no VEV 
but large Yukawa couplings yr to the leptons/neutrinos, 
while Hu2 could have large VEV but small Yukawa Yj^ . 
The contribution from Fig. 1 is given by 

4 

{mN'=)ab ^ (ig^2-)3 "^i/2 - I0^^^mi/2tan^ P , (17) 

which is sub-eV as long as the model lies in the low 
tan /? = K'2/ Ki regime. 

B. Collider Signatures 

The virtue of the LR model, is the possible lepton num- 
ber violation at the colliders through the production of 
Wr . This happen when the RH neutrino picks up a large 
Majorana mass, and the end result is same-sign dilep- 
tons [l7j . This generic feature of non-supersymmetric 



LR theories becomes somewhat subtle here. First, the 
heavy RH neutrino A'^i is tied to the extra neutralino Z' 
as given Eq. ([4]), where Mz' is now given by Eq. (7). 
Therefore, if the Dirac mass dominates, the RH neutrino 
Ni is not necessarily lighter than the Wr gauge boson 
as assumed in the conventional case, and the dilepton 
signal would be somewhat suppressed. Second, even if 
A'^i is made lighter by taking a fairly large toi/2 > Mz', 

the mixing between Ni and Z' , equal to Mz' /mi/2 can- 
not be too small, for the sake of naturalness. In this 
case, one generally expects A^i to decay through a slep- 
ton which could in principle be on shell (Fig. [3]) . Possible 
multi-lepton signatures from such processes are discussed 
separately in Ref [isj . 




FIG. 2: Production and decay of the RH neutrino. 



There is a potential problem regarding the possibly low 
scale left-right symmetry in this model. The point is that 
the scalar fields with no VEVs lead to tree-level flavor- 
changing processes, and must be very heavy, on the order 
of lOTeV [1^. There is a possibility of fine-tuning these 
effects against the usual supersymmetric ones [l^] , but it 
goes against the principle of low-energy supersymmetry 
being behind naturalness. The trouble is that the mass 
of the heavy Higgs is obtained from the D-term, and thus 
roughly equal to the Wr mass. The most natural out- 
come is then to have all the new scales raised to 10 TeV, 
which would take them out of the LHC reach. Of course, 
some of the soft mass terms could be made smaller as 
to bring some superpartner masses to the TeV region. 
In any case, barring cancellation, it is unlikely to have 
Wr^ accessible at the LHC in this model. This would 
not change any of our conclusions regarding the neutrino 
masses and our point of two light neutrinos would remain 
equally valid, even for very heavy Wr. 
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C. Comparison with other LR models 

Before concluding this section, we briefly comment on 
the supersymmetric LR model with the additional Higgs 
fields. It obviously makes no sense to introduce more 
doublets on top of the slcptons, which favors the model 
with triplet Higgs. In this case, the RH scale is not au- 
tomatically connected with the scale of supersymmetry 
breaking. One can use higher dimensional operators to 
break SU{2)r at some large scale (the AMS model) plf . 
The interesting feature of this theory is the necessarily 
exact R-parity and the stability of the LSP [l^] . Further- 
more, one expects doubly charged supcr-multiplcts to re- 
main light enough to be accessible at the LHC [21, 23 2^ . 

On the other hand, if the higher dimensional opera- 
tors were absent or very small in this minimal model 
with Higgs triplets, one would be forced to break R- 
parity (the KM model) [2^ . Otherwise, one ends up 
breaking electromagnetic charge invariance. Just as in 
the minimal model studied here, the RH sneutrino must 
also develop a non-vanishing VEV, in addition to the 
triplet Higgs VEV. The difference between the FS and 
KM models lies in the RH neutrino spectrum. In the 
KM model (as in the AMS one), all the RH neutrinos 
are made heavy by the triplet Higgs VEVs, so it might 
be less appealing in view of the recent BBN and neutrino 
oscillation progress. As a payback, the KM model pos- 
sesses richer collider phenomenology. In the KM model, 
there is another contribution to the multi-lepton signal 
through 7V2.3 and the e'^ — mixing [l^. Therefore, 
when Ni is made heavier than Wn (if mi/2 ^ "^zOj such 
signal will be turned off in the FS model while it could 
still happen in the KM model. 

IV. SUMMARY 

B — L symmetry plays a special role in the SM: it is 
the only anomaly free global symmetry. This strongly 
suggests its gauging, which requires the existence of RH 
neutrinos, one per each family. In turn, one predicts 
non- vanishing neutrino masses through the seesaw mech- 
anism. In the MSSM, B — L symmetry can play another 
important role in setting the stage for the R-parity. Re- 
cently, important progress has been made in the con- 
struction of the supersymmetric B — L model, by getting 
rid of any additional Higgs beyond the MSSM ones. As 
a consequence, the B — L scale is tied to the supersym- 



metry breaking scale, and if one believes that low energy 
supersymmetry is accessible at the LHC, so is the B — L 
gauge boson. 

In this letter, we find another important characteristic 
feature of this theory: besides the usual actives neutri- 
nos, one predicts automatically two light sterile neutri- 
nos with masses in the sub-eV region. This is the central 
result of our paper. It is equally valid in the left-right 
extension of the theory, where B — L gauge symmetry 
follows automatically. We further discussed some new 
possible phenomenological features of the LR theory rel- 
evant for LHC. 

It was recently argued these new light states are fa- 
vored by the BBN. The crucial point is that the B — L 
gauge interactions assure them being abundant at the 
time of nucleosynthesis. They also may be playing an 
important role in neutrino oscillations, as long as there 
is an oasis in the desert above the B — L scale, around 
10^ GeV or so. 
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Note added After this paper was submitted to the arXiv, 
we learned from Rabi Mohapatra that the idea of using 
RH sneutrino VEV in order to break the SU (2)^ symme- 
try is originally due to [2^ , and in the case of B — L it was 
used in (23 |. Furthermore, it was noticed in j29j that only 
one RH neutrino gets a large Majorana mass which is the 
starting point of the analysis presented here. This point 
is carefully discussed in [30| , which appeared on the same 
day as our paper. We wish to acknowledge too that the 
possibility of neutrinos having both Majorana and Dirac 
mass terms on competing level was discussed recently in 
the context of the supersymmetric B ~ L theory in [3l| . 
However, the neutrino spectrum and phenomenology is 
different from the one in this work. 
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